The 6j-symbol: Recursion, Correlations and Asymptotics 
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We study the asymptotic expansion of the {6j}-symbol using the Schulten-Gordon recursion rela- 
tions. We focus on the particular case of the isosceles tetrahedron and we provide explicit formulas 
for up to the third order corrections beyond the leading order. Moreover, in the framework of spin- 
foam models for 3d quantum gravity, we show how these recursion relations can be used to derive 
Ward-Takahashi-like identities between the expectation values of graviton-like spinfoam correlations. 

o . 

O 1 I- INTRODUCTION: {6J}-SYMBOL AND THE RECURSION FORMULA 

m 

Spinfoam models present a proposal for a well-defined path integral formulation for quantum gravity. They define 
the structure quantum space-time at the Planck scale. The current challenge is to study their semi-classical behavior 
' at large scale, show that we recover general relativity and a locally flat space-time as expected and then to extract 
rj 4 - the perturbative quantum gravity correction to the classical gravitational dynamics. The main spinfoam models for 
4d quantum gravity are the Barrett-Crane model [l[ and the more recent family of models @, H[ based on coherent 
state techniques. These models are all derived from the reformulation of general relativity as a constrained topological 
BF-theory and attempt to define a discretization of the path integral over space-time geometries. 

A recent proposal, which is actively investigated, to probe the semi-classical behavior of the spinfoam amplitudes 
is the test of the graviton propagatorproposed by Rovelli and collaborators [f|. There have been several more or 
less rigorous analytical computations , as well as numerical simulations 0] , showing that the most recent spinfoam 
models have the expected behavior reproducing at first order (and for the simplest space-time triangulation) Newton's 
law and the tensorial structure of the graviton. These rely on the computation of the leading order asymptotics 
(N ' of some spinfoam amplitudes for large spin. The next step is to go beyond the leading order and compute the 
quantum corrections which should correspond to the loop corrections of the standard perturbative quantum field 
t-H | theory approach. This necessarily requires a better understanding of the structure of the asymptotics of the spinfoam 
0^ . amplitudes. Here, we investigate this issue in the simplified framework of 3d quantum gravity. 

In three space-time dimensions, gravity is a topological theory which can be exactly quantized as a spinfoam 
model, the well-known Ponzano-Regge model !8|. Its basic building block is the {6j}-symbol. Its asymptotics has 
been well-studied at leading order and been derived using various techniques: from the brute-force calculation based 
on the explicit formula of the {6j}-symbol in term of factorials [9( to more refined saddle point approximation on 
integral formulas [Io| . Recently, these calculations have been pushed further in order to compute the corrections to 
the asymptotic behavior using both the saddle point technique [ll[ and the brute-force method [l2|. These results 
apply to the computation of graviton-like correlations in 3d quantum gravity [ll], [l4| . This should lead to a better 
understanding of the structure of the quantum corrections of the spinfoam graviton propagator, which should be 
relevant to the four-dimensional case. 

In the present paper, we are interested in the computation of the asymptotics of the {6j}-symbol through the use 
of recursion relations. As was first shown in the {6j}-symbol satisfies a recursion formula which are intimately 
related to the topological invariance of the Ponzano-Regge spinfoam model. This formula turned out to be very useful: 
it can be approximated in the large spin limit by a second order differential equation and one can use it to derive the 
leading order of the asymptotics through a WKB approximation, but it also allows fast numerical calculations of the 
{6j}-symbol. 

Here, we investigate two aspects of these recursion relations. First, we show how to extract the next-to-leading 
and subsequent corrections to the {6j}-symbol. We focus on the case of the isosceles tetrahedron (since this is the 
case relevant for the computation of graviton-like correlations [3) and we compare our results with the previous 
work [ll|, [13] • Second, we study the consequences of the existence of such a recursion relation on the behavior of 
the graviton-like correlation functions in spinfoam models. We show that it leads to relations on these correlations 
functions, which relate the expectation values of different observables. These relations are similar to the Ward- 
Takahashi identities (and to the Schwinger-Dyson equation) in standard quantum field theory, which allow to relate 
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different correlation functions of different orders. These equations play a key role in the study of the renormalization 
of quantum field theories, so we expect these new recursion relation to be equally relevant to the study of the 
renormalization/coarse-graining of spinfoam models. 



II. THE RECURSION RELATION FOR THE ISOSCELES TETRAHEDRON 



A. Exact and Approximate Recursion 

We will focus on the isosceles tetrahedron, which is relevant for the computations of geometrical correlations in the 
simplest non-trivial toy model in 3d quantum gravity [14| . Such a tetrahedron has four of its edges of equal length 
with the two remaining opposite edges of arbitrary length. The corresponding isosceles {6j}-symbol is: 



{a,b}j = 



a J J 
b J J 



where J £ N/2 and a, b are integers smaller than 2 J (to satisfy the triangular inequality). The associated tetrahedron 
has edge lengths lj = dj/2 for j — a, b, J, where dj = 2j + 1 is the dimension of the SU(2) representation of spin j. 
The volume V of the tetrahedron is given by the simple formula: 



Vj(a,b) = -lj b ^4lj-(ll + q), 



(1) 



while the (exterior) dihedral angles 9j can also be written in term of the edge lengths (see e.g. [Tl( for more details): 
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The general recursion relation for the {6j}-symbol given by Schulten and Gordon in |13j simplifies in this specific 
isosceles case: 



(la- 



4lj 



(la 



1 



{a + !,&}„ 



21, 



(4l 2 T -l 2 a )cos(e a ) + 



1. 



{a,b}j + (l a --) 4lj-(l a --) 2 {a-l,b}j = Q 

In the asymptotic regime, we know (analytically and numerically) the behavior of the {6j}-symbol at the leading 
order: 



{a,b}j ~ {a, b} L j 



LO — 



1 



Vl2nV 



(laOa 



+ l b 6 b + Al.,0., + - 



(4) 



which is actually valid under the assumption that the tetrahedron with edge lengths l a ,h,lj exists (else the generic 
asymptotics can be expressed in term of Airy functions). The oscillatory phase is given by the Regge action Sr — 
l a Q a + hdb + 4lj6j. Using the obvious trigonometric identity cos((n + 1)4>) + cos((n — l)</>) = 2cos0 cosrw^, we can 
write an exact recursion relation for the leading order of the {6j}-symbol: 



y/Vj(a + l,b){a + 1, b} L j° - 2 cos 9 a y/Vj(a,b){a, b} L } ° + y/Vj(a - l,b){a - 1, b} L 7 



LO 



0. 



(5) 



A similar recursion relation holds for ^-shifts and also J-shifts. 

The most natural idea is to compare this recursion relation for the leading order to the previous equation on the 
exact {6j}-symbol to see how to use them to extract the next-to-leading correction to the asymptotic behavior. We 
can first find the link between the leading order of equation Q and the leading order of equation ([5]). Both equations 
can be written under the same form at the leading order: 



{a+ l,b}j - 2cos9 a {a,b}j + {a - !,&}„ 



0. 



(6) 



which turns into a simple second order differential equation in the large spin limit. Then the next-to-leading order of 
the equation ((3]): 



y/Vj(a, b) (l + ^ (l - apjj)) {a + 1, b}j - 2 cos 6 a ^Vj(a, b){a, b}j 

/TTl TT 1 /-, 2l 2 „ x\ r 



+v^M) 1- Ha-l,6}j»0 



(7) 



will have to be compared to an recursion relation for the next-to-leading order of the {6j}-symbol. 
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B. Pushing to the Next-to-Leading Order 

We are interested in the asymptotic expansion of the {6j}-symbol. It was shown in previous works [lj], 03 that lj 
seems to be the right parameter to consider when studying the semi-classical behavior of the {6j}-symbol. So from 
now we write: 



a J J 
b J J 



= {la, h}lj- 



Notice that shifting a by ±1 is equivalent to shifting the edge length Z a = <x+l/2by±l. We rescale now lj by XI j 
and we replace the exact {6j}-symbol by a series in 1/A alternating cosines and sinus of the Regge action (shifted by 
7r/4) in the previous equation © . The fact that there is no mixing up of cosines and sinus at all order was show in 
[111 ]. More precisely, we write the {6j}-symbol asymptotic expansion under the form: 

{Xl a , Xl b } Mj = x3/2D ^ lM) [cos(XS R + tt/4)+ f(1>( '- Wj) sin(AS fl + tt/4) + Gll Hia.h,h) cos(XS b . + tt/4)) 



.E^*M C0S (XS R 



•tt/4) 
tt/4)- 
-tt/4) 



tt/4)) 



Gl3) ";/ t ' ij) cos(g fl + V4) 

G(4) V"' j) sin (^ + ^/ 4 ) + °( A ~ 5 )]> 



(8) 

where the pre-factor denominator D(l ai l , lj) is given by the square-root of the tetrahedron volume as in equation ([4]). 
To study the asymptotics, it is convenient to factorize the whole equation ^ by A 3 / 2 . We then write {l a ± 1/A,4}ij 
for {Xl a ± 1, Xli,}xij. We also factorize the coefficients of the recursion relation. We start by defining C(lj) = 

'a(4/j - l 2 a ) = 16(A( ' a ^ J '' j))2 where A(a,b,c) = \yj{a + b + c)(a + b - c)(a - b + c)(-a + b + c) is the area of the 
triangle of edge lengths given by a, b and c. The coefficient which appears in front of {l a ± 1/A, lb]ij becomes 
C{l a ± 1/(2A), l b , lj) = {l a ± 1/(2A))(4/ 2 T - (l a ± 1/(2A)), where we underline that the shift is ±1/(2A) and not simply 
±1/A. We expand C(l a ± 1/(2A), lb, lj) in term of derivatives: 



1 1 r) n C 

^±1/(2A),M,) = £^^ W 



with 



c = i a (M*j-il) 



3li 



d-c 

d 3 C 
91" 



= -U 2 a 

= -6 

for n > 4 



(9) 



Then to express {l a ± 1/X,l b }ij we need to expand D(l a ± 1/A), F^{l a ± 1/A), and G (l) (^ ± 1/A): (i 6 {!• ••4}) 



D{l a ±l/X) = D±\I£ 
^)(/ Q ±l/A) = Et l o(-l) 
GW(/ a ±l/A) = E 4 fe r (-l) 



1 d 2 £> 
\k 1 d fc F<'> 



_i_ 1 d 3 D 

31 sD^'air 



fc!A fc azj 
\fc 1 d fc G w 
fe!A fe dl* 



1 8 4 £> 



(10) 



F^(lj) was computed in a previous paper [HI, [13]. It was also suggested that the asymptotic expansion of the {6j}- 
symbol in term of the length scale A is given by an alternative of cosines and sinus at each order, so we expect that 
G^\lj) = for Vi > 1. Finally, we also need to expand the Regge action XSji(l a ± j), remembering that 9j = Oj(l a ) '■ 



XSfi(l a ± — ) — XSr 



« (_l)fc+i d k 6a 
2^(k + l)\X* dl* 



(11) 
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with 

dOg _ -2lgl b 

dl a (■iPj-lD^iPj-ll-ll 

d 2 9 a _ 2l b (4l 2 l 2 a -2l 4 a -l 2 J 2 b + 16l i J -4l 2 b l 2 J ) 
~dTf ~ (4'j-i^) 2 [4i'j-/2-i^]3/2 

d 3 9 a _ 

~~di*~ — (u'j-ii) 3 yu 2 ,-ii~ii} 5/2 

= W-ia*n&,-i*-i>)v* M 8l a + 8l « l b + IWjWb ~ ™ljl a + 7l 6 Jt + 35201% - U72ljl a l* + 21% 
+14QZ*#» - 560^5 - 3840/«/ 2 + 2432/6/2/2 + ^a^a + 2048 ^2 _ 448/ 6 /4 _ 3072Z io + 32 #4)f 6 ) 

We can now write an asymptotic recursion equation from equations ©, ©, ©, (fTT))) and (fTTj) in terms of A neglecting 
terms of order 0(A -4 ) and smaller, assuming that A is large. This leads to a couple of equations at each order, one 
for the cos-oscillations and one for the term in sin: 

• The first equation is given by the terms of order A and it is trivially satisfied (0 = 0) since we have already written 
the leading order of the {6j}-symbol proportional to cos(Sr + j) (the Ponzano-Regge asymptotic formulae). 



• The second equation is given by the terms of order A 



-l. 



1 dC 1 dD\ . ,„ , 100. 



sm 



2Cdl a Ddl a J y ' 2 dl a 
which can be rewritten as a differential equation for D: 



--^cos(0 a )=O (12) 



dlnD 1 



dl a 2 



d6 a cos 9 a d In C 



dl a sin a dl a 



(13) 



This allows to determine D: \nD = | ln(C sin(0 o )) + K, which simplifies into D = K\l l a h^^lj — l% — l£ where 
if is a constant factor. Thus this second equation shows that D is correctly proportional to the square-root of 
the volume V of the isosceles tetrahedron. To determine the normalization constant K (as well as G^), the 
orthonormality property of {6j}-coefficients can be employed: ^ a Al a \/hlb' {a, b}j{a, b'}j — Sw and we get the 
K = \J\2-k. The details are given in the next section. 



The third equation is given by the terms of order A 2 and which are proportional to cos(5r + j) 

d-F (1) _ la i ( J_dC_ ]_dD_\ l§6a_ i 1 d 2 9 a 

\2C dl a DdlaJ 2 dl a 6 dll 



AC sin 8 a 

I COS 9 a 



( 1 d 2 D 3_9fC _i_J_£D/'j_aC_J_£D_i_ /l*^)^ 

\2D~ST a r 8C ~7H^ D dla \ 2C dl a D dl a ' \2 dl a > )) 



(14) 



where we used the fact that ^^fp — 7)§r) sm (^a) + lfl s " cos (^a) = ^ (eqn. ITS"!) ) to remove all the terms 

proportional to itself. The first term of the right-hand side of the equation (|14p comes from the variation 
of the coefficient in front of {a,b}j in the recursion equation ([3]). The terms with a derivative of C with respect 
to l a come from the coefficients in front of {a ± 1, b} j and {a, b}j. The variation of C with respect to l a is given 
by the variation of the areas of the triangles of the tetrahedron. From eqn. (fl"3|) . we relate it to the variations 
of D (the volume) and to the variations of the dihedral angle 9 a : Tjfp = ^fp — sing" W L ' ^ ne terms with a 
derivative of D with respect to l a come from the variation of the leading order of the asymptotic of the {6j}- 
symbol and the terms with a derivative of the dihedral angle 9 a come from the variations of the Regge action 
Sr. We can now compute the derivative of with respect of l a (equation (TT4|) ) in terms of l a , h and lj the 
edge lengths of the tetrahedron: 



^WT = ~ 48(^(-^ + ^) 2 (^-^-^)( S /^ ib ) (- 32 ^j^ + 10#* + 96/^ - 960/6/4 + 15 ^4 + 400| 4j4 f a _ 100/^^4 
-l&l 2 / a ll - 1664/f/ 2 / 2 + 20/8/ 2 + 576/j/ 2 /4 + 3072/*/ 2 - 3072/}° + A8ljl a + 2304/«/ 2 - 5761%) 

(15) 

and then easily integrate this equation over / a : 
= 

24(^^)(2^^)(4^^^) 3/2 / a i 6 +*Vb,tJ) 
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The integration constant Z(lb, Ij) can be determined using the symmetry properties of the {6j}-symbol: sym- 
metry of the isosceles {6j}-symbol with respect to l a and Z&, coupling of l a , lb and lj by this isosceles {6j}-symbol 
and homogeneity of ([F^] = lj 1 ) imply that Z(lb,lj) — 0. Then this gives us the same result as in the 
previous paper [ijj ■ Moreover, using the definitions of the tetrahedron volume ([T]) and of the dihedral angles 
(J2J), we can express F^ in terms of some geometrical characteristics of the tetrahedron: 



48(12^)3/^8 1 j 

The fourth equation is given by the terms of order A -2 and which are proportional to sm(Sn + f )■ It is the 
same equation as the previous one for G^ 1 -* but the right-hand side is now equal to zero (homogenous equation). 
That is we simply get that 

so G^ 1 -* = Z(lt,,lj) is just a constant of integration. Once again the symmetry properties of the {6j}-symbol 
implies that G (1) = 0. 

The next equation is given by the terms of order A -3 and which are proportional to sin(S*fl + § )■ We get an 
equation for the first derivative of F^ 2 ' with respect to l a 

dF^(ij) _ c OS e a a 2 ^ 1 ' _ ( i ajg , of' 1 ' 

dl a 2 sin a dl'i V sin 2 B a dl a ' J dl a 

( -LOR J_ §Q. , 1 d 2 P _ ( J_dD_\ 2 1 d 2 C , J_ / 1 de^"] 2 \ 1 d6„. , ( J_dP_ 1 

I D dl a 2C dl a 2P dl 2 \P dl a ) 8C dl 2 a + 3! ^ 2 dl a J J 2 dl a ^ \D dl a 2C 

f J_ dC_ _ J_dP_\ 1 ( ld6a\ 2 i 1 d 2 9g 1 09 a , / 1 3c\ 
^ 2C Oia D dl a J 2 \2 dl a J ^ 3! "STf - 2 91, \ ^ <9i a / 



^ (3(12^)8 - (nv)Xit (m - i b ) 2 + Kif) - «) + ®l 2 i 12 



sin 6 a 



1 dC 1 d 2 D I 1 8 2 C 1 8J) i / 1 8C 181)11(1 99 a \ , 1 d 2 g a 1 39 a , / 1 dD \ 1 dD 1 d 2 D , l_a 3 D 

2C 3i„ 2D "511" t 8C "SlJ" D Oio V 2C dl a D dla I 2 \ 2 dla I + 3! 2 dl a \ D dla I D dla D ~WT 3<D ' 

* 2 



■3\C dl 3 n 



J_dC_ ( 1 dD\ 

2C dla \P dla ) 



(19) 

We recall that D is proportional to the square root of the tetrahedron volume, G can be expressed in terms of 
the volume V and the sinus of the dihedral angle 6 a (see equation (TT3l)V To integrate this equation, we first 
express explicitly Ht in terms of l a , Zf, and Zj, and then deduce F^: 

F {2 \h) = 4608 ( (4 ^^^) 2 ( 4 ^^g 1 )2(4 p j „ z g^ ; g ) 3^ f g ) (-2359296ZgZ^Zg - 224512^,4 + 100/ l2j4 + 57^4,12 + 112896 ^| 

+2727936^/ 1 / + 5308416^ 6 + 212$°*^ - 5898240Z 2 J* 4 - 11520^ o i^ + 941056^Z^6 + 31584 0j Z 6 
-2416Z^io/2 _ 79872/8/^2 _ 480Z 12 /2/2 _ 7040/8/ 6 / 2 , _ 24m 1 a ljtf + 100l 4 Jl 2 + 212/1 / 6 , 
+2727936/^ - 700416/^/g - 5898240/^ 4 / 2 - 11520/^° - 700416/^1° + 609/8/8 + n 2 896/8/8 
+576/1 2 / 4 - 2359296/ 4 /i°/ 2 + 528384/^/8/2 + 5849088/ 2 /} 2 / 2 - 79872/ 2 / b 8 / 6 7 + 31584/ 4 /8/5 - 7040/6/8/2 7 
-224512/4/6/6 + 58816/6/6/4 + 8640/i°/5/ 2 + 8640/ b 10 /4/ 2 - 480/i 2 / 2 / 2 + 528384/ 2 /8/ 6 ) 

which is the only result with the required symmetries 2 . The geometrical meaning of this function does not seem 
obvious. Nevertheless, we can give a more compact expression for the denominator of F^ 2 ': 

(4/ 2 - Z 2 ) 2 (4Z 2 - Z 2 ) 2 (4Z 2 - Z 2 - llfllll = (21) 

COS uj 



1 ^(h) = - 2 3 04((4^- 1 g) 1 3 ( 4 lS - l ;)3 (4iaj - 1 g +l g)4 1 g) (- 1604f g'g'j + 1250816/4/^6 _ 207104^6/f/e + 3Wm u>t*li - 16934«J[igI« - 
3920^i 10 ^ + 24992i6«8« 4 T - 46848iJ°i^? - 7129088i4« 1 / °i 4 + 1770496i«Z|« 4 + 16832«^i°« 4 T + 34368«|«5 J 6 ~ 6816 ^ 2 ' " 278912i«Z^4 + 
14524416(^ 1 r 2 ( 4 + 486144«2i8;8 - 560iJ 2 i 4 (2 _ 43776i 2 Zi°Z5 - 3317760Za 2 ZJ 10 Zg + 22241280i 4 (i 2 i 2 + 794Z 42 Zf - 6955008/^ W + 
2801664«i 2 Z| + 672iJ 4 i 2 i 2 - 26542080« 4 i^ 4 - 451584^°^ + 46080«^J° - 2304i 12 ^ _ 21233664? 1 / + 37158912« 2 ii 6 + 1072128Z|i^b ~ 
1528ia 12 / 4 i 2 - 10911744^ 4 « 4 - 35979264i 2 ;i 4 / 2 + 1728«l 2 « 4 .« 2 + 9953280^i 1 / 2 + 228096«i°«| - 10368Z^ - 2073600«|i 1 ,0 + 27Z 4 °Z° + 
23592960(/z 2 + 400/^1° - 88ll 4 lf - 814411° iffi + 1001™ 1%) 

2 If the result is not symmetric after integration, a non-null integration constant has to be added and its determination can be done using the 

symmetry properties of the {6j}-symbol. Indeed, we have — = H(l a ,l b ,lj) so by integration over l a , F^(lj) = h(l a ,l b ,lj)+Z(l b ,lj). 
Moreover by symmetry, we must have g t = H(l a = = la,lj) and then integrating over we obtain a second expression for 

p( 2 ) : F^{lj) = h(l a = lb,h — la,lj) + Z{la,lj) which implies that the constant of integration satisfies Z(l b ,lj) — Z(l a ,lj) = h(l a = 
Zf,,/;, = l a ,lj) — h(l a ,l b ,lj). This equation allows to determine Z and to get (1206 . 
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• The next equation comes from the terms of order A 3 which are proportional to cos(>!?r + j): 

-OL^'O < 22 » 

which implies once again that G^ 1 = Z(l b , Ij) is a constant of integration. Then the symmetry properties of 
the {6j}-symbol implies G [2 \l,) = 0. 

We can now give the asymptotic expansion of an isosceles {6j}-symbol until the next to next to leading order 
(NNLO): 

{L, «r LO = . =| = : = \cos(S R + F« { l 3 ) sin (S R + + F^ fe) cos(S R + j)} (23) 

where the expression for F^ and F^ are given by equations (fT7|) and (|2"0)) . This result seems to confirm that the 
expansion of the {6j}-symbol is a series alternating cosines and sinus of the Regge action (shift by j). In the case of 

an equilateral tetrahedron, all the edges have the same length, that is l a = h = l.J = I and V = -jjl 3 - Then equation 
([2"3"| reduces to: 

{6j}^° = ; 1 cos(S fl + -) == sm(S R + -)- 1 —= cos(S R + -) (24) 

q VW2 4 722!/4 2 5/2V^ 1 4 ; 20736 2 V*2V^ V 4 ; 17 

where the Regge action is given by S R = 610 and 8 = 9 a = 9 b = 0j = arccos( — 1/3). This result is confirmed by 
numerical simulations. The plot figure [T] represents numerical simulations of the equilateral {6j}-symbol minus its 
approximation Moreover, to enhance the comparison, we have multiplied by t 1 1 2 to see how the coefficient of the 
NNLO is approached and we have divided by sin^^ + j) (oscillations of the next to next to next to leading order) 
to suppress the oscillations. This gives an error that decreases as expected as l^ 1 . 

0.006 ' 
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0.004 

0.003 

0.002 

0.001 ' . 



1000 2000 3000 4000 5000 



FIG. 1: Difference between the equilateral {6j}-symbol and the analytical result (|24p . The x-axis stands for I = d/2 and d goes 
from 100 to 5000. The error decreases as expected as Z -1 confirming our asymptotic formula. 



• The next two equations come from the terms of order A 4 . The equation for G^ is the same as the one for G^ 
and G^ 2 \ that is: - = 0. Using the same arguments of symmetry we deduce that G^ — 0. This confirms 



our expectation of a series alternating cosines and sines in the asymptotic of the {6j}-symbol: 



{Xl a , Mb}\ij — 



X 3 / 2 D{l a ,l b ,lj) 



cos{XS R + -) 



2^ Tfc cos(XS R + - + e(k)~) 

fc=i 



(25) 



where e(k) = — 1 when k is odd and e(fc) = when k is even. We already have an expression for F^ and F^- 2 \ 
The equation for F^ is the second equation of order A -4 and gives its first derivative with respect to l a in 
terms of l a , l b and lj. It is straightforward (though lengthy) to integrate it over l a and get the expression 3 of 
F^ in terms of l a , l b and lj. In the equilateral case (l a =h = lj = I), the formula reduces to: 



i?(3) 



28833535 



1 



17915904 2 9 / 2 2V4p 



(26) 
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Therefore, we have the expression of the asymptotic expansion of the equilateral {6j}-symbol up to the next- 
to-next-to-next to leading order (NNNLO): 



{6j} 



NNNLO 



cos(S R + -) - if^rjjj sm(S R + -) 

45673 7T, 28833535 . ti\ 

- 2073624?2 cos(5« + -) + 1791590429/2/3 ™P* + 4 , 



(27) 



We check this result numerically by computing it using Mathematica for two values of spins j = 50 and 

({67} -{6i} NNNLO )Z 9/2 

j = 100. More precisely, we computed the renormalized error ° q " c ' os (s R +^) an< ^ we S ^ the expected 

1/A-behavior. However for I > 100, Mathematica is not accurate enough and the numerical errors are too 
important to get exploitable results. In the general isosceles case, the expression of is quite complicated 
and its geometrical interpretation remains to be understood. Nevertheless, we can again give as before a more 
compact formula for the denominator of : 

denominator^, = 3317760(4*5 - l 2 a - l 2 b f/ 2 l 3 jUtf ~ llf^j ~ = 30(48) 3 (28) 

From this equation, the equation giving the denominator of F^ and remembering that the denominator of F^ 
can be written under a similar form: denominator^) = 48 cos ^ g , we can conjecture that: 

. (12Vj(a,b)) 3k 
denominator^ <fc) cx — , — (29) 



(cos ( 



where F^ are the terms appearing in the asymptotic expansion of the {6j}-symbol (|25|) . And consequently, 
the numerator of is a polynomial in lj of degree 8k. 

So, using the recursion relation for the isosceles {6j}-symbol as well as its symmetry properties, we have computed 
explicitly the asymptotic expansion of the isosceles {6j}-symbol to the fourth order up to an overall factor K (this 
integration constant K comes from the integration of the first equation (fT3|) ). The well-known value K — \/Ylix which 
already appears in the Ponzano-Regge formula can be obtained easily using the unitary property of the {6j}-symbol, 
as we show in the next section. The equilateral case has been checked against numerical calculations. This method 
using the recursion relation is fairly easy to implement. It requires integrating a rational fraction at each level and 
does not involve neither Riemann sum nor saddle point analysis. Moreover, since the coefficient C of the recursion 
relation © is a polynomial of degree 3; = for n > 4. Therefore, we expect to get a stable relation for the first 

derivative of F^ and G^ k > with respect to l a for k > 3. On one hand, this allows to prove that G^ always vanishes; 
and on the other hand, it should provide a systematic method to extract F^ for arbitrary order k. 

We conclude this section with a general remark on the asymptotic expansion of the {6j}-symbol. In the context 
of 3d quantum gravity, it is often argued that the leading order of the {6j}-symbol is a cos(S) instead of a complex 
phase exp(+iS), thus reflecting that the path integral is invariant under a change of (local) orientation (see e.g.[15j). 
This obviously neglects the +7r/4 shifts, which can be considered as a quantum effect (like an ordering ambiguity). 
However, in the light of the present expansion, it is clear that we have terms of the type sin(5) beyond the leading 
order and such terms are not invariant under the change S — ► —S. This means that the role of this symmetry in the 
spinfoam path integral should be more subtle than originally thought. 



- 28270080ZJ 6 1% 



1 F^th) = r-s-, 3 k, , — . — ,,, , ^ , — =-^(-6965741568/?! + 1069728/^ ifl 2 , 

3317760(ll(l a -2lj)(l a +2lj)(4,l 2 J -ll) 3 (4,l 2 J -ll-ll)(S/2)ll(4l 2 J -i 2 .) 2 ) v a J b ' a b J 

7430401$% ll 8 + 1750503260160Z 2 / 3 Z 2 ll + 379404800zj 2 fil? + 788705280/1° if I 4 - 33547200^i 4 Z 4 -98578608Zi 2 Z*Z5-81032970240«^ 4 Z 
379404800«^ 2 Z 6 7 - 389214720Z* Z 2 ^ 4 - 323813376000^ - 156036/^J 6 - 22176Z 2 Z 4 Zj s _ 3824640Z* 6 Zj,Z 4 - 28408^ ' K 
1262270545920« 1 r s «?iJ - 1262270545920Z 1 T 8 Z 4 i 2 . + 461814497280Z 1 T 6 Z 2 Z? + 1069728ii 6 2 2 ,Z 6 . + 74144841728/5 l\ 2 lt + 8 



.„ if -31104000* - 
8245202780162 1 /« 4 Zf + 
6861127680W|/^ 6 - 
+509607936000Z 24 - 
' 2 Zj 2 - 28408M 8 



7430401™ Ijl? - 22176ll s l 2 jl 4 - 81032970240Z| / 4 « 2 zf + 7201368Z, 

1401753600« i ; 2 «\ 2 -267541217280«| l r 4 i^g+568565760i^5'b + 111 
3981312000^° If - 156036«J 6 if -511602^° Zi 4 +344217600Zj 4 if - 14 
3981312000Z^ 4 Zj° - 511602Zi 4 « 1 ° - 31104000(J 6 Z| - 9852420096002 
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C. Consequences of the unitary property of the {6j}-symbols 



The orthogonality property of the {6j}-symbols states that: 

^ 4^aV Ibh' {lg, h}lj{lg, h'}lj = ^l b l b , 



(30) 



This relation corresponds to the unitarity of the evolution in the Ponzano-Rcgge 3d quantum gravity. We want to 
use this property to determine the constant of the leading order of the {6j}-symbol. From the recursion relation we 
have shown that {l a , h}f® — i K cos(Sr + j); but K is still undetermined. For large spin and for lb 4', we 

can approximate the unitary property at the leading order in (lb — lb>) by: 



Jo 



K 2 7T 7T 

dlattJb T/ , cos(S R (l a , l b ) + -) cos(S R (l a , l b ') + ~r)~ S(l b - l b >). 
Vj(a,b) 4 4 



(31) 



The product of the cosines can be simplified at leading order: 



cos(S R (l ai l b ) + -) cos(S R (l a , l b 



1- 
4' 



1 

2 
1 

2 

il 2 j-ll-2l 2 a 



cos(S R (l a ,l b ) + S R (l a , l b ') + - ) + cos(S R (l a , l b ) - S R (l a , l b ')) 
cos(2S R (l a , h) + |) + cos((Z fc - h')6 b ) 



where the dihedral angle 6b = arccos ^— Al -' / J ! h _ l t a J is considered as a function of the length l a . We do a saddle point 

approximation. The first term oscillates and its integral is exponentially suppressed. And we are left with the second 
term, which should satisfy the following equation: 

K 2 



j — : 



dlalalb'-, 



We recall that: 



therefore we can conclude that 



Vj(a,b) 

i r 00 

— / dl a cos(l a (l 
2?r J_oo 



cos((l b - h>)0b) ~ S(l b - lb') 



(32) 



b - h')) = S(lb - h') 



b V ~ 27T 



d0 b 



dla 



(33) 



9b and l b are so conjugate variables and K comes from the Jacobian of the change of variables between l a and 6b- 
Computing the derivative of the dihedral angle gives: 

36b -2 -lj b „ 1 

dl a 



6Vj(a,b) 



K = 



T2tt 



(34) 



Moreover, pushing the approximation of the unitary property to the next to leading order in (4 — lb> ) and using the 
next to leading order of the {6j}-symbol shows that = 0. This was already shown in the previous part using the 
recursion relation and the symmetry properties of the {6j}-symbol and comes as a confirmation. 



III. "WARD-TAKAHASHI IDENTITIES" FOR THE SPINFOAM GRAVITON PROPAGATOR 

We are interested in the two-point function in 3d quantum gravity for the simplest triangulation given by a single 
tetrahedron. This provides the first order of the "spinfoam graviton propagator" in 3d quantum gravity. 
Considering the isosceles tetrahedron, we focus on the correlations between the two representations a and b: 

(0(a)0(b))^ = | J2 Mo)Mb)0(a)0(b){a, b} J} Z = ^ ^j(a)^j(b){a, b}j, (35) 

a, b a, b 

where ipj(j) is the boundary state, which depends also on the bulk length scale J, and O, O are the observables whose 
correlation we are studying. 

Now, inserting a recursion relation with shifts on a, b or J in the sum over the representation labels b leads to 
equations relating the expectation values of different observables. We distinguish two cases: when the state -0J does 
not change or when the length scale J also varies. 
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A. Relating Observables 

Inserting the recursion relation on a-shifts in the definition of the correlation function, we obtain the following exact 
identity: 



(^f^O(a - l)0(b)(l a - ±)(4Z 2 - (/„ - I) 2 ))^ -(0(a)0(&)2Z Q (2cos(? a (4Z 2 - Z 2 ) + i))^ 

H^f^Oia + l)OWa + U^J (la + \?))^ = 



We call this a Ward identity for our spinfoam correlation. If the observable diverges at a = 0, more precisely if 
it contains terms in 1 ja or in l/(a + 1), then we need to take into account extra boundary terms in this equation 
corresponding to contributions at a = 0. But all observables usually considered are regular in this sense. 

Then one can choose different sets of observables O and O and one gets different identities on the correlation 
functions of the spinfoam model. For example, taking 0(a) — l a , we get: 

(^fer^W^ " - V2)(4Z5 - (l a - l/2) 2 ))„ -(d(&)(2cos^ 2 (4Z 2 - l\) + P a /2))^ 

+< 3 *$5?0(&)(l + l)(l„ + l/2)(Alj - (l a + 1/2) 2 )),, = 0. 



We recall that the area of the triangle of edge lengths given by l a , lj, lj is equal to A(l ai lj) = jl a ^j41 2 j — Z 2 ; then 

(/„ ± l)(l a ± l/2)(4/ 2 - (/„ ± 1/2) 2 ) = 16[A 2 (l a ± 1/2, lj) ± ^n^UA ], therefore we can rewrite the previous 

equation as an equation between correlation functions of the observable 0(b) and different observables proportional 
to the square of the triangle area A(l a , lj): 

(Mf^l[A 2 (l a - 1/2, lj) - A J ( Ci/i j) P(b)U -{{2cose a A\l a ,lj) + ll/2)d{b))^ 

+( t ^S 1 {^(i a + 1/2, ij) + A j^iJi 3) w)h = °- 

The standard choice of boundary is a phased Gaussian 0, 0, [l6[ : 

MJ) ~ e^e- 2 ^^, (37) 

where $ is a fixed angle defining a posteriori the external curvature of the boundary and a is an arbitrary real positive 
number (which can be fixed by the requirement of a physical state 16]). In this case, we can compute explicitly the 
ratios ip(a ± l)/ip(a) entering the Ward identity: 

ipj(a± 1) ■ - - --- '"-'■> 



Of course, this ratios does not depend on 6; therefore if the observable 0(b) = 1, then the dependence on b only 
appears in one correlation function through the cosine of the dihedral angle 9 a . As another example, we consider 

0(a) = la 1 and 0(b) = then: 

/ ^j(a-l) l a -l/2 4l 2 j-(l a -l/2) 2 4l 2 j-ll \ 9 / rriq fl U-l-t I 1 U 2 ,-ll \ 

\ ^Jj(a) l a -l IP, 475 z \ Luf,t7 a^7 4/7 r J^T 4 ,a /V 

I / ^j(o+l) iq + 1/2 (4i 2 7 -(i a + l/2) 2 4#-ig v _ n 

which can be approximated by: 

(e -*2 V ^ A(a a -l/2) 2 )Aa 2 ))^-2 e ^(cos0 a A(/ 2 )Aa 2 ) + ^A(Z 2 ))^ + ( e ^ e - 4tt ^ i A((Z a +l/2) 2 )A(Z 2 ))„ « 



where A(i|) = 



;:-4ij 



B. Rescaling the Tetrahedron 

We can now vary also the length scale lj. First let's notice that in the same way we wrote an exact recursion 
relation for the leading order of the isosceles {6j}-symbol shifting the representation a (equation (JSJl), we can write a 
similar exact recursion relation for the leading order of the {6j}-symbol shifting the label </; that is 

VK/ +1 (a,6){a, 6}^ - 2 cos(40j) ^Vj(a, b){a, b} L j° + A /y J _ 1 (a,h){a, &}£i = (38) 
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Inserting this recursion relation on J— shifts in the definition correlation function, we obtain the following identity: 

-2(cos(40j) VK/(a, b)0(a)0{b))^ =0 (39) 

The correlation functions appearing in this equation are in fact approximation. We are allowed to use the leading 
order of the {6j}-symbol because the boundary state used picks the function on large jo- And for the same reason, 
we can expand y/V J±1 (a,b) and the ratios ^^$)tp^\b) '■ 



(l - ^l\i_ vi ) e 4Q J ~"° " [1+ "A"j-°--* l .> ] 0(a)0(b))i, - 2(cos(40j)^VA^b)O(a)O(b))^ 

r — / n; \ -In ( 2l J-('a + 'b)) M 31 j -2(l a + i b ) , „ 

+ < v A0M) (l + 4 ^!' ; f_^ ) e 4Q ~, [1 ~^.J-^ ] 0{a)0{b)) 4 , » 0. 



(40) 



We hope that such equation will turn out useful to study the asymptotic properties of the correlations function as 
the length scale J grows large, but we leave this for future investigation. 



Conclusion 



We have used the recursion relation satisfied by the {6j}-symbol to study the structure of its asymptotical expansion 
for large spins. The exact recursion relation allowed us to compute explicit the asymptotical approximation of the 
isosceles {6j}-symbol up to fourth order. This confirms previous results [111, |T3] and introduces techniques allowing 
further systematic analytical calculations of the corrections to the behavior of the{6j}-symbol at large spins. However 
a clear and simple geometrical interpretation of the polynomials appearing in this expansion is still missing, but the 
differential equations that we provide for these coefficients should be a first step in this direction. 

This work is useful in particular for the study of large scale correlations in the spinfoam model for 3d quantum 
gravity. In this context, the recursion relation allowed us to write equations satisfied by the spinfoam correlations 
similar to the Ward identities of standard quantum field theory. We hope that such recursion techniques can be 
further applied to the study of 4d spinfoam amplitudes and the resulting spinfoam graviton propagator [17[ ■ 
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